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Abstract 

We show that some common and important global 
constraints like All-Different and GCC can be 
decomposed into simple arithmetic constraints on 
which we achieve bound or range consistency, and 
in some cases even greater pruning. These de- 
compositions can be easily added to new solvers. 
They also provide other constraints with access to 
the state of the propagator by sharing of variables. 
Such sharing can be used to improve propagation 
between constraints. We report experiments with 
our decomposition in a pseudo-Boolean solver. 



1 Introduction 

Global constraints allow users to specify patterns that 
commonly occur in problems. One of the oldest and most 
useful is the All-Different constraint |Lauriere, 1978[ 
Regin, 19941. This ensures that a set of variables are 
pairwise different. Global constraints can often be de- 
composed into more primitive constraints. For example, 
the All-Different constraint can be decomposed into 
a clique of binary inequalities. However, such decompo- 
sitions usually do not provide a global view and are thus 
not able to achieve levels of local consistency, such as 
bound and domain consistency. Considerable effort has 
therefore been invested in developing efficient propagation 
algorithms to reason globally about such constraints. For 
instance, several different propagation algorithms have been 



developed for t he All-Di ffere nt constraint l|Regin, 1994 
Leconte, 19961 |Puget, 1998| [Mehlhorn and Thiel, 2000" 
Lopez-Ortiz et al, 2003|. Inthis paper, we show that several 



important global constraints including All-Different can 
be decomposed into simple arithmetic constraints whilst 
still providing a global view since bound consistency can be 
achieved. 

There are many reasons why such decompositions are in- 
teresting. First, it is very surprising that complex propagation 
algorithms can be simulated by simple decompositions. In 



many cases, we show that reasoning with the decompositions 
is of similar complexity to existing monolithic propagation al- 
gorithms. Second, these decompositions can be easily added 
to a new solver. For example, we report experiments here us- 
ing these decompositions in a state of the art pseudo-Boolean 
solver. We could just as easily use them in an ILP solver. 
Third, introduced variables in these decompositions give ac- 
cess to the state of the propagator Sharing of such variables 
between decompositions can increase propagation. Fourth, 
these decomposition provide a fresh perspective to propagat- 
ing global constraints that may be useful. For instance, our 
decompositions of the All-Different constraint suggest 
learning nogoods based on small Hall intervals. 

2 Formal Background 

A constraint satisfaction problem (CSP) consists of a set of 
variables, each with a finite domain of values, and a set of 
constraints specifying allowed combinations of values for 
some subset of variables. We use capitals for variables and 
lower case for values. We write dom{X) for the domain 
of possible values for X, min{X) for the smallest value 
in dom{X), max{X) for the greatest, and range(X) for 
the interval [min{X),max{X)]. A global constraint is one 
in which the number of variables n is a parameter. For 
instance, the global All-Different([Xi, . . . , con- 
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straint ensures that Xi ^ Xj for any i < j I Regin, 199^) . 
We will assume values range over 1 to d. 

Constraint solvers typically use backtracking search to ex- 
plore the space of partial assignments. After each assignment, 
propagation algorithms prune the search space by enforc- 
ing local consistency properties like domain or bound con- 
sistency. A constraint is domain consistent (DC) iff when a 
variable is assigned any of the values in its domain, there exist 
compatible values in the domains of all the other variables of 
the constraint. Such an assignment is called a support. A con- 
straint is bound consistent iBC) iff when a variable is assigned 
the minimum or maximum value in its domain, there exist 
compatible values between the minimum and maximum do- 
main value for all the other variables. Such an assignment is 
called a bound support. Finally, between domain and bound 
consistency is range consistency. A constraint is range con- 
sistent (RC) iff when a variable is assigned any value in its 
domain, there exists a bound support. 

Constraint solvers usually enforce local consistency after 



each assignment down any branch in the search tree. For this 
reason, it is meaningful to compute the total amortised cost 
of enforcing a local consistency down an entire branch of the 
search tree so as to capture the incremental cost of propaga- 
tion. We will compute complexities in this way. 

3 All-Different constraint 

The All-Different constraint is one of the most useful 
global constraints available to the constraint programmer For 
instance, it can be used to specify that activities sharing the 
same resource take place at different times. A central concept 
in propagating the All-Different constraint is the notion 
of a Hall interval. This is an interval of m domain values 
which completely contains the domains of m variables, [a, b] 
is a Hall interval iff \{i \ dom{X,) C [a, b]}\ = b ~ a + 1. 
In any bound support, the variables whose domains are con- 
tained within the Hall interval consume all the values in the 
Hall interval, whilst any other variables must find their sup- 
port outside the Hall interval. 

Example 1. Consider an All-Different constraint over 
the following variables and values: 



Xi 

X2 
X3 

Xi 

X5 

[1, 1] is a Hall interval of size 1 as the domain of 1 vari- 
able, X5 is completely contained within it. Therefore we can 
remove [1, l]from the domains of all the other variables. This 
leaves X2 with a domain containing values 2,3,4. 

[3, 4] is a Hall interval of size 2 as it completely contains 
the domains of 2 variables, Xi and X^. We can thus remove 
[3, A\from the domains of X2 and X4. This leaves the follow- 
ing range consistent domains: 

1 2 3 4 5 



Xi 
X2 
X3 
Xi 
X5 

Enforcing bound consistency on the same problem does not 
create holes in domains. That is, it would leave X4 with the 
values 2,3,4,5. 

To identify and prune such Hall intervals from the domains 
of other variables, Leconte has proposed a RC propagator for 
the All-Different constraint [Leco nte, 1996| that runs in 
Q{n^) time. We now propose a simple decomposition of the 
All-Different constraint which permits us to enforce RC. 
The decomposition ensures that no interval can contain more 
variables than its size. We introduce 0{nd^) new 0/1 vari- 
ables, Aiiu to represent whether Xi takes a value in the inter- 
val [/, u]. For l<i<n, l</<u<(i and u — I < n, we 
post the following constraints: 



E 



< 



X., G [l,u] 
u-l + 1 



(1) 
(2) 



Example 2. Consider again the last example (i.e. an 
All-Different constraint on X I e [3,4], X2 G [1,4], 
X3 e [3, 4], X4 e [2, 5] and X5 e [1, 1];. 

First take the interval [1, 1]. Since X^ G [1, 1], (|7) implies 
A511 = 1. Now from (|2]l, X]i=i — 1- That is, at most 
one variable can take a value within this interval. This means 
that A211 — 0. Using ([7} and A211 = 0, we get X2 ^ [1, 1]. 
Since X2 G [1, 4], this leaves X2 G [2, 4]. 

Now take the interval [3,4]. From (|7}, A134 — ^334 = 1. 

Now from (|2|, X)i=i ^i34 ^ 2. That is, at most 2 vari- 
ables can take a value within this interval. This means that 
A234 = A434 = 0. Using (0 we get X2 ^ [3, 4], ^ [3, 4]. 
Since X2 G [2,4] and X4 G [2,5], this leaves X2 = 2 
and Xi G {2, 5}. Local reasoning about the decomposi- 
tion has thus made the original All-Different constraint 
range consistent. 

We will prove that enforcing DC on the decomposition en- 
forces RC on the original All-Different constraint. We 
find it surprising that a simple decomposition like this can 
simulate a complex propagation algorithm like Leconte's. In 
addition, the overall complexity of reasoning with the decom- 
position is similar to Leconte's propagator 

Theorem 1. Enforcing DC on constraints (|7} and (O en- 
forces RC on the corresponding All-Different constraint 
in 0{nd^) down any branch of the search tree. 

Proof: jLeconte, 1996| provides a necessary and sufficient 
condition for RC of the All-Different constraint: every 
Hall interval should be removed from the domain of variables 
whose domains are not fully contained within that Hall inter- 
val. Let [a, b] be a Hall interval. That is, \H\ = b — a + 1 
where i7 = {i \ dom{Xi) C [a, 6]}. Constraint ([T]! fixes 
Aiab — 1 for all i G i?. The inequality with I = a and 
u — b becomes tight fixing Aiab = for all i ^ H. Con- 
straint ^ for I ^ a, u — b, and i ^ H removes the interval 
[a, b] from the domain of Xi as required for RC. 

There are 0{nd^) constraints ^ that can be woken 0{d) 
times down the branch of the search tree. Each propagation 
requires 0(1) time. Constraints ([T]i therefore take 0{nd^) 
down the branch of the search tree to propagate. There are 
0{d^) constraints (|2|i that each take 0{n) time to propagate 
down the branch of the search tree for a total of 0{nd^) 
time. The total running time is given by 0{nd^) + 0{nd^) = 
0{nd^). □ 

Note that if we use a solver in which we can specify that 
constraints only wake on reaching a particular bound, we can 
decrease to 0(1) the number of times constraints ([T]i are wo- 
ken, which gives a total complexity in 0{nd^). 

What about bound consistency of the All-Different 
constraint? By using a representation that can only prune 
bounds [Ohrimenko et ai, 2007 ], we can give a decomposi- 
tion that achieves BC in a similar way. In addition, we can 
reduce the overall complexity in the case that constraints are 
woken whenever their bounds change. We introduce new 0/1 
variables, Bik, I < k < d and replace ([T]l by the following 
constraints: 



We illustrate this decomposition on our running example. 



Bu = 1 
Aiiu — 1 



X, < I (3) 
- A - 1) (4) 



Theorem 2. Enforcing BC on constraints Q fo en- 
forces BC on the corresponding All-Different constraint 
in 0{ncP) down any branch of the search tree. 
Proof: We first observe that BC is equivalent to DC on 
constraints (|2|i because Au^ are Boolean variables. So, the 
proof follows that for Theorem[T] except that fixing Au^ = 
prunes the bounds of dom{Xi) if and only if = 

or Biu = 1, that is, if and only if exactly one bound of the 
domain of Xi intersects the interval [l,u]. Only the bounds 
that do not have a bound support are shrunk. The complexity 
reduces as (O appears 0{nd) times and is woken 0{d) times, 
whilst (01 appears 0{nd^ ) times and is woken just 0(1) time. 
□ 

A special case of All-Different is Permutation 
when we have the same number of values as variables, and 
the values are ordered consecutively. A decomposition of 
Permutation just needs to replace Q with the following 
equality where (as before) 1 < I < u < d, and u — I < n: 



E- 

i=l 



1 



(5) 



This can increase propagation. In some cases, DC on con- 
straints ([TJ and (|5]) will prune values that a RC propagator for 
Permutation would miss. 

Example 3. Consider a PERMUTATION constraint over the 
following variables and values: 

12 3 



X2 

Xs 

These domains are range consistent. However, take the inter- 
val [2,2]. By DC on (Q, A122 = ^222 = 0. Now, from 0, 
we have ^i22 = 1- Thus A322 — 1- By ([T]), this sets 

X3 = 2. On this particular problem instance, DC on con- 
straints ((7) and (O has enforced domain consistency on the 
original All-Different constraint. 

4 GCC constraint 

A generalization of the All-Different con- 
straint is the global cardinality constraint, 

GCC([Xi, . . . , X„], [h, . . [Ui,.. .,Um]). This 

ensures that the value i occurs between Zj and Ui times in Xi 
to Xn. The GCC constraint is useful in resource allocation 
problems where values represent resources. For instance, in 
the car sequencing problem (probOOl at CSPLib.org), we can 
post a GCC constraint to ensure that the correct number of 
cars of each type is put on the assembly line. 

We can decompose GCC in a similar way to 
All-Different but with an additional O(d^) integer 
variables, Niu to represent the number of variables using val- 
ues in each interval [I, u]. Clearly, Niu E ^i' J2^=i 
and Nid = n. We then post the following constraints for 
l<i<n,l<l<u<d,l<k<u: 



A, 



. = 1 ^ 
Niu = 



Xi e [l,u] 



1=1 



'A 



ilu 



Nik + N^k+i)u 



(6) 
(7) 
(8) 



Example 4. Consider a GCC constraint with the following 
variables and upper and lower bounds on the occurrences of 
values: 



V 


1 


2 


3 


4 


5 


Xi 


* 










X2 




* 


* 






Xs 












Xi 




* 








X5 














1 


1 





1 


1 




5 


5 


5 


5 


5 



Enforcing RC removes 1 and 3 from X2,X4 and X^ and 
leaves the other domains unchanged. We can derive this from 
our decomposition. From the lower and upper bounds on the 
number of occurrences of the values, we have Na £ [1,5] 
except for N^'^ £ and we have N12 G [2,10],A^i3 G 

[2, 15] and Nu e [3, 20]. By A333 = 1. From ([71, 



33 = Y.i=i ^i33 e [1, 5]. From N15 = Nu + N55 we have 
A^i4 G [3,4] (i.e., upper bound decreased) because Ni^ = 5 
and N 55 G [1,5]. Similarly, we derive from Nu = ^'^13+^^44 
that Ni3 G [2, 3] and from N13 = iVi2 + A'^33 that N12 G 
[2, 2]. /toot the same constraint, we shrink N13 to [3, 3] and 
N33 to [1, 1]. Finally N12 = A^n + A^22 shrinks Nn to [1, 1]. 
By (0, Am = ^1333 = 1, so by 0, A.n = 0, i, G 2. .5 and 
^i33 = 0, i, G {1, 2, 4, 5}. By (|6l), this removes 1 and 3 from 
X2, X/i, X^. Local reasoning about the decomposition has 
made the original GCC constraint range consistent. 

We next show that enforcing DC on constraint (|6]l and BC 
on constraints O and ([8]l enforces RC on the GCC con- 
straint. 

Theorem 3. Enforcing DC on constraint ^ and BC on con- 
straints and dS)) achieves RC on the corresponding GCC 
constraint in 0(nd^) time down any branch of the search tree. 

Proof: We use ly for the number of variables Xi whose 
range range{Xi) intersects the set V of values, and Sy for 
the number of variables Xi whose range is a subset of V. We 
first show that if RC fails on the GCC, DC on Q and BC on 
O and ([8]l will fail. We derive from iQuimper et al., 2005 



Lemmas 1 and 2] that RC fails on a GCC if and only if 
there exists a set of values V such that 5*^ > J^vev '^^ 
such that Iv < Tliv^v^""- Suppose first a set V such that 
Sy > J2vev'^v The fact that domains are considered as 
intervals implies that either range{V) includes more vari- 
able domains than the sum of the upper bounds (like V), 
or the union of the range{Xi) that are included in V lets 
a hole of unused values in V, which implies that there ex- 
ists an interval [/,u] C V such that „] > Tlveli u]'^v 
So, in any case, there exists an interval [/, u] in V with 
S[i,u] > EvelLu]^^'- By ® we have ^"^1 A,i„ > S^^^^^ 
whereas the greatest value in the domain of Niu was set to 
T,ve[i,u] ■ So BC will fail on Niu = J2"=i ^Uu- Sup- 
pose now that a set F — {vi, . . . ,Vk} is such that ly < 
J2v ev ^vi - The total number of values taken by Xi variables 
being equal to n, the number of variables Xi with range{Xi) 
not intersecting V is greater than n — J^v ev^vi, that is. 



Thanks to (Q, we know that for any /, u, Niu > S^^ 

Nl(vi-l) +iV(t,i + l)(t,2-l) ^ 



„i- So, 



+ N{vk+i)d > n-Y^^.^ylv 



The initial domains of Niu variables also tell us that for 
every Vi in V, Ny.y. > ly.. Thus, mm(-/Vi(„^_i)) + 
TOm(iV„^„J+mm(A^(„j+i)(„2„i)) + . . .+mm(A^(„^+i)d) > 
n = max{Nid)- Successively applying BC on Niy-^ — 
7Vi(.„i_i)+A^«it,i,thenon7Vi(„3_i) = iVi„, +Ar(„^+i)(^„_i), 
and so on until A^i^ = iVi^^. + N(^Vk+i)d will successively 
increase the minimum of these variables and will lead to a 
failure on A'^i^. 

We now show that when DC on Q and BC on ^ and (O 
do not fail, it prunes all values that are pruned when enforcing 
RC on the GCC constraint. Consider a value v G dom{Xq) 
for some g e l..n such that v does not have any bound sup- 
port. We derive from | Quimper et al. , 2005 Lemmas 1 and 



6] that a value v for a variable Xq does not have a bound sup- 
port on GCC if and only if there exists a set V of values such 
that either (i) Sy — X^wiey v &V and range{Xq) is not 
included in V, or (ii) ly — J^wev lw,v^V and range{Xg) 
intersects V. In case (i), V contains v and the values it con- 
tains will be taken by too many variables if Xq is in it. In case 
(ii), V does not contain v and its values will be taken by not 
enough variables if Xq is not in it. Consider case (i): Since 
DC did not fail on (|6]l, by a similar reasoning as above for de- 
tecting failure, we derive that V is composed of intervals [I, u] 



such that S, 



[Lu] 



Consider the interval 



containing v. The greatest value in the initial domain of Ni 



lu 



was 



we[Lu] ^w- 



which is exactly the number of variables 



with range included in [I , u] without counting Xq because its 
range is not included in V . Thus, O forces Aqi^ — and (|6]l 
prunes value v from dom{Xq) because v G [/, u] by assump- 
tion. Consider now case (ii): V = {wi, . . . , Wfc} is such that 
ly = X^tj ey ■ The total number of values taken by the Xi 
variables being equal to n, the number of variables Xi with 
range{Xi) not intersecting V is equal to n — J2v £V ^vi, that 

-1] + '^[t.i + l,i;2 -1] + • ■ • + ^[vk + i ,d] = - Y.v,<£V ^-"i ■ 

Thanks to (|7]i, we know that for any I, u, Niu > ^j. So, 

^i(«i-i)+iV(t,i+i)(„2-i)+. . ■ + > n-J2v,ev ^vi- 

The initial domains of iV;„ variables also tell us that for 
every Vi in V, N^.y. > ly.. Thus, mm(A^i(i,j_i)) + 
min{Ny^yJ+min{N(y^+i)(y^_i)) + . . .+min{N(y^+i)d) > 
n — max{Nid). Successively applying BC on Niy^ = 
^^i(i,i-i)+^^t.ii)i,thenon7Vi(^2_i) = A^it,i +iV(„i+i)(i,2-i), 
and so on until Niy^ = A^^-^j,.!-) + Ny^y^ will increase 
all min(A^i(„._i)) and min{Niy.), to the sum of the min- 
imum values of the variables in the right side of each con- 
straint so that TOin(iVii,^,) = min{Ni{^y^_i))+min{Ny^y^) + 
min(iV(„j+i)(„2-i)) + . . . + min{Ny^y^). Then, because 
max{Nid) = n, BC on Nid = iVi^^ + iV(„j^+i)d will de- 
crease the maximum value of Niy^ and iVj-^j.+i)^ to their 



minimum value, BC on Niy^ — N- 



+ Ny 



will de- 



crease the maximum value of Ni(^y^_i) and Ny^y,^ to their 
minimum value, and so on until all -/V(t,;+i)(i,._^i_i) are forced 



to the singleton 



At this point, (Q forces A 



= s, 



[■Ui + 1,D, + 1-1]- 

for every 



variable Xj with range not included in the interval 



in S, 



1] because that interval is saturated by variables Xp 
for which = 1. By as- 



sumption value V is not in V, so there exists such an interval 
[vi + l,Vi+i — 1] that contains v. Furthermore, range{Xq) 
intersects V, so it is not included in [vi + 1, — 1]. There- 
fore, A 



is forced to and (|6]l prunes v from 



q(vi + l)(vi+i-l) 

dom(Xq). 

There are 0{nd^) constraints ^ that can be woken 0{d) 
times down the branch of the search tree in 0(1), so a total 
of 0{nd^) down the branch. There are 0{d^) constraints (Q 
which can be woken 0{n) times each down the branch for 
a total cost in 0{n) time down the branch. Thus a total of 
0{nd^). There are 0{d'^) constraints ([8]l that can be woken 
0{n) times down the branch. Each propagation takes 0(1) 
time to execute for a total of 0{nd^) time down the branch. 
The final complexity down the branch of the search tree is 
therefore 0{nd^) + 0{nd^) + 0{nd^) = 0{nd^). □ 

What about bound consistency of the GCC constraint? As 
in the case of All-Different, by replacing constraints (|6]l 
by constraints (O and (|4]i, the decomposition achieves BC. 

Theorem 4. Enforcing BC on constraints 0, ^ and 
dSj achieves BC on the corresponding GCC constraint in 
0{nd^) time down any branch of the search tree. 

Proof: The proof follows that for Theorem [3] except that fix- 
ing Aiiy — prunes the bounds of dom{Xi) if and only if 
exactly one bound of the domain of Xi intersects the interval 
[l,u]. The complexity reduces to 0{nd^) as BC on (|3]l and 
(IDi is in 0{nd^) (see Theorem|2]i and BC on ^ and ^ is in 
0{nd'^) (see Theorem[3]l. □ 

The best known algorithm for BC on GCC runs in 0{n) 
time at each call | Quimper et al. , 2005| and can be awaken 
0{nd) times down a branch. This gives a total of O^n^d), 
which is greater than the 0{nd^) here when n > d. Our 
decomposition is also interesting because, as we show in the 
next section, we can use it to combine together propagators. 

5 Other global constraints 

Many other global constraints that count variables or values 
can be decomposed in a similar way. For example, the global 
constraint Same([Xi, . . . , X„], [Yi, . . . , y„]) is satisfied if 
and only if the Yi variables are a permutation of the Xi vari- 
ables. A monolithic flow-based propagator for this constraint 
is given in jBeldiceanu et al., 2004| . The following decom- 
position encodes the Same constraint where 1 < i < n, 
l<l<u<d, l<k and k < u: 



Aiiu — 1 



X, G 



Bily, — 1 



G 



E 

i=l 



Nly 



Ik 



obtained by posting de- 
.,X„],[Oi,...,0„i]) and 
) and eliminating com- 



This decomposition can be 
compositions for eGCC ( [Xi , . . 

EGCc([yi,...,y„],[Oi,...,o„ 

mon sub-expressions (eGCC is an extended form of the GCC 
constraint in which upper and lower bounds on occurrences of 
values are replaced by integer variables). This is another ar- 
gument in favor of decompositions since it allows constraints 
to share "internal" state through common intermediate vari- 
ables. Such sharing can increase propagation. 



Example 5. Consider the following example: 



2 3 4 5 



Y2 



If we have Oi G [0, 1] for 1 < « < 5 then 
both eGCC([Xi,X2], [01,02,03,04,05]) and 
EGCC([ri,r2], [01,02,03,04,05]) are BC. 

However, enforcing BC on the decomposition of 
Same([Xi, X2], [11, ^2]) removes 3 from the domain 
of X2 and Yi- 

In fact, we conjecture that enforcing BC on this decom- 
position achieves BC on the Same constraint itself. Similar 
decompositions can be given for other global constraints like 
NValue and Common. 



6 Experimental Results 

To test these decompositions, we ran experiments on 
pseudo-Boolean encodings (PB) of CSPs containing 
All-Different and Permutation constraints. We used 
the MiniSatH- 1.13 solver on an Intel Xeon 4 CPU, 2.0 Ghz, 
4G RAM with a timeout of 600 seconds for each experiment. 
Our decompositions contain two types of constraints: Sum 
constraints like ^ and MEMBER constraints like ([TJ. The 
Sum constraints is posted directly to the MiniSatH- solver. 
To encode MEMBER constraints, we use literals Bij for the 
truth of Xi < j lOhrimenko et ah, 2007[ , and clauses of the 
form {Aiiu — \) (i3i(;_i) = A 3^ = 1). This achieves 
bound consistency (Theorem|2]i. To increase propagation, we 
use a direct encoding with literals Zij for the truth of Xi — j 
and clauses {Auu = 0) {Zij = 0), j € [I, u]. The overall 
consistency achieved is therefore between BC and RC. We 
denote this encoding HI. To explore the impact of small 
Hall intervals, we also tried HIk, a PB encoding with only 
those constraints (|2]i for which u — I + 1 < k. This detects 
Hall intervals of size at most k. Finally, we decomposed 
All-Different into a clique of binary inequalities, and 
used a direct encoding to convert this into SAT (denoted BI). 



Pigeon Hole Problems. Table [T| gives results on pigeon 
hole problems (PHP) with n pigeons and n — 1 holes. Our 
decomposition is both faster and gives a smaller search tree 
compared to the BI decomposition. On such problems, de- 
tecting large Hall intervals is essential. 

Double-Wheel Graceful Graphs. The second 
set of experiments uses double-wheel graceful 
graphs l |Petrie and Smith, 20031. We converted the CSP 
model in | |Petrie and Smith, 2003[ into a PB formula. This 
model has an All-Different constraint on node labels 
and a PERMUTATION constraint on edge labels. For the 
Permutation constraint we use (|5]l. We strengthen the 
BI decomposition with clauses to ensure that every value 
appears at least once. Table |2] show that our decomposition 
outperforms the augmented BI decomposition on many 
instances. Whilst detecting large HaU intervals can greatly 



n 


BI 
bt/t 


Hh 
bt/t 


HI3 
bt/t 


bt/t 


bt/t 


HIg 

bt/t 


5 
7 
9 
11 
13 
15 
17 
19 
21 


30/ 0.0 
622/ 0.0 
16/35/ 0.3 
998927/ 29.3 

-/- 
-/- 
-/- 
-/- 
-/- 


28/ 0.0 
539/ 0.0 
18455/ 0. / 
665586/44.8 

-/- 
-/- 
-/- 
-/- 
-/- 


4/ 0.0 

47/ 0.0 
522/ U.U 
5681/0.3 
13876/ 0.9 
1744765/ 188.6 
-/- 
-/- 
-/- 


6/ 0.0 

1 Tl / A A 

122/ 0.0 
171/ 0.0 
2568/ 0.2 
24109/2.6 
293762/ 48.0 
107780/21.8 
-/- 


8/ 0.0 

180/ 0.0 
247/ 0.1 
1054/ 0.2 
8989/ 1.1 
857175/368.0 
550312/ 426.2 


10/0.1 
195/ 0.1 
165/ 0.1 
4219/ 0.6 
39713/ 9.9 
57817/ 33.5 



Table 1: PHP problems, t is time and bt is the number of 
backtracks to solve the problem. 



BI 
bt/t 



176/0.1 
30/ 0.1 
22/ 0.2 

1341/ 1.0 
2948/ 3.6 
2418/5.5 
3378/ 8.6 
19372/ 118.3 
839/ 5.4 



bt/t 



90/0.1 
14/0.1 

526/ 0.4 
873/ 0.9 
2047/ 4.2 
724/ 2.2 
1666/5.7 
9355/66.: 
12356/84.2 



HI3 
bt/t 



214 



63/ 0.1 

212/0.2 
87/ 0.3 
318/ 0.7 

1710/3.6 
643/ 2.8 
1616/9.0 
120/85.9 
1556/ 13.9 



Hh 
bt/t 



1290/ 1.7 
1212/2.9 
1574/4.0 
368/ 2.4 
30/ 1.8 
10/ 2.1 
14/ 2.4 



Hh 
bt/t 



27/ 0.9 

3955/ 19.5 
10123/ 129.7 
4051/35.0 
7456/ 105, 



Hh 
bt/t 



405/ 6.5 
5709/71.2 
25552/92.7 



Table 2: Double-wheel graceful graphs, t is time and bt is the 
number of backtracks to solve the problem 



reduce search, in some cases the branching heuristics ap- 
pear to be fooled by the extra variables introduced in the 
encodings. 

Overall these experiments suggest that detecting Hall inter- 
vals reduces search significantly, and focusing on small Hall 
intervals may be best except on problems where large Hall 
intervals occur frequently. 

7 Other Related Work 

The All-Different constraint first appeared in the AL- 
ICE constraint programming language (Lauriere, 1978[ . Re- 
gin proposed a DC propagator that runs in 0{n''-^) time 



i Regin, 1994| . Leconte gave a RC propagator based on Hall 
intervals that runs in O(n^) time |Leconte, 19961. Puget 
then developed a BC propagator also based on Hall intervals 
that runs in 0(nlog(n)) time I Puget, 1998|. This was later 
improved by Melhorn and Thiel IjMehlhorn and Thiel, 2000| 
and then Lopez-Ortiz et al. | Lopez-Ortiz et al., 2003) . 

The global cardinality constraint, GCC was introduced 
in the CHARME language l Oplobedu ef a/., 1989| . Re- 
gin proposed a DC propagator based on network flow 
that runs in O(n^) time iRegin , 1996| . Katriel and 
Thiel proposed a BC propagator for the eGCC con- 
straint iKatriel and Thiel, 2 003 1 . Quimper et al. proved 
that enforcing DC on the eGCC constraint is NP-hard 
I Quimper ef a/., 20041. They also improved the time com- 



plexity to enforce DC and gave the first propagator for en- 
forcing RC on GCC. 

Many decompositions have been given for a wide 
range of global constraint. However, decomposition 
in general tends to hinder propagation. For instance, 
I Stergiou and Walsh, 19991 shows that the decomposition of 



All-Different constraints into binary inequalities hinders 



propagation. On the other hand, there are global constraints 
where decompositions have been given that do not hinder 
propagation. For example, Beldiceanu et al. identify condi- 
tions under which global constraints specified as automata 
can be decomposed into signature and transition constraints 
without hindering propagation | |Beldiceanu et al, 2005) . 
As a second example, many global constraints can be de- 
composed using Roots and Range which can themselves 
often be propagated effectively using simple decom- 
positions [Bessiere et al, 2005[ [Bessiere et al, 2006at 
[Bessiere et al, 2006b) . As a third example, decompositions 
of the Regular and Cfg constraints have been given 
that do not hinder propagation iQuimper and Walsh, 2006; 
Quimper andW alsh, 200 7] [Quimper and Walsh, 200^ 
Bessiere et al, 2008, ,Katsirelos et al, 2008| . As a fourth 
example, decompositions of the Sequence constraint have 
been shown to be effective l |Brand etal., 20071 . Finally, the 
Precedence constraint can be decomposed into ternary 
constraints without hindering propagation | |Walsh, 2006| . 

8 Conclusions 

We have shown that some common global constraints like 
All-Different and GCC can be decomposed into simple 
arithmetic constraints whilst still maintaining a global view 
that achieves range or bound consistency. These decompo- 
sitions are interesting for a number of reasons. First, we 
can easily incorporate them into other solvers. Second, the 
decompositions provide other constraints with access to the 
state of the propagator. Third, these decompositions provide 
a fresh perspective on propagation of global constraints. For 
instance, our results suggest that it may pay to focus propa- 
gation and nogood learning on small Hall intervals. Finally, 
these decompositions raise an important question. Are there 
propagation algorithms that cannot be efficiently simulated 
using decompositions? In [Bessiere etal, 2009], we use cir- 
cuit complexity to argue that a domain consistency propaga- 
tor for the All-Different constraint cannot be simulated 
using a polynomial sized decomposition. 
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